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Abstract 

We present the results of an experimental investigation of the acoustics and 
fluid dynamics of Tibetan singing bowls. Their acoustic behavior is rationalized 
in terms of the related dynamics of standing bells and wine glasses. Striking 
or rubbing a fluid-filled bowl excites wall vibrations, and concomitant waves at 
the fiuid surface. Acoustic excitation of the bowl's natural vibrational modes 
allows for a controlled study in which the evolution of the surface waves with 
increasing forcing amphtude is detailed. Particular attention is given to ratio- 
nalizing the observed criteria for the onset of edge-induced Faraday waves and 
droplet generation via surface fracture. Our study indicates that drops may be 
levitated on the fiuid surface, induced to bounce on or skip across the vibrating 
fiuid surface. 

1 Introduction 

Tibetan singing bowls are thought to have originated from Himalayan fire cults 
of the 5th century BC and have since been used in various religious ceremonies, 
including shamanic journeying and meditation. The Tibetan singing bowl (see 
Fig. [T]) is a type of standing bell played by striking or rubbing its rim with a 
wooden or leather- wrapped mallet. This excitation causes the sides and rim 
of the bowl to vibrate and produces a rich sound. Tibetan bowls are hand 
made and their precise composition is unknown, but generally they are made of 
a bronze alloy that can include copper, tin, zinc, iron, silver, gold and nickel. 
When the bowl is filled with water, excitation can cause ripples on the water 
surface. More vigorous forcing generates progressively more complex surface 



wave patterns and ultimately the creation of droplets via wave breaking. We 
here quantify this evolution, and demonstrate the means by which the Tibetan 
singing bowl can levitate droplets. 



Figure 1: Our Tibetan singing bowls: a) Tibet 1, b) Tibet 2, c) Tibet 3 and d) 
Tibet 4. 

Related phenomena are known to appear in other vessels, including the 
Chinese singing bowl, on which vibrations are generated by rubbing the vessel's 
handles with moistened hands. The more familiar vibration of a wine glass is 
produced by rubbing its rim with a moist finger. The dependence of the glass's 
vibration frequency on its material properties, geometry, and characteristics of 
the contained fiuid was elucidated by French J] and subsequent investigators [H 
|3l m [5] . The coupling between two singing wine glasses has been investigated by 
Arane et al. [6] . Apfel [7^ demonstrated experimentally that wine glass vibration 
generates capillary waves near the walls of a fluid-filled glass. Moreover, he made 
the connection between these waves and the acoustic whispering gallery modes 
elucidated by Rayleigh ^ . By studying the deformation and the sound spectrum 
produced by a single wine glass, Rossing [9 elucidated the mechanism of the 
glass harmonica [10], an instrument designed by Benjamin Franklin. Joubert et 
al. pT provided a theoretical rationale for observations of standing waves in a 
singing wine glass. The Tibetan singing bowl has to date received relatively little 
attention. Inacio et al. examined experimentally the acoustic response of bowls 
excited by impact and rubbing [12 , and proposed a dynamical formulation of 
the bowl and presented some numerical simulations [13 . The hydrodynamics 
of a fluid-filled Tibetan bowl will be the focus of our investigation. 

In 1831, Faraday [14j demonstrated that when a horizontal fiuid layer is 
vibrated vertically, its interface remains flat until a critical acceleration is ex- 
ceeded. Above this threshold, a fleld of waves appears at the interface, para- 
metric standing waves oscillating with half the forcing frequency [151 IISl Ull HI 
[T9l |2Q] . The form of such Faraday waves depends on the container geometry; 
however boundary effects can be minimized by using a large container. The 
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Faraday waves have a wavelength prescribed by the relative importance of sur- 
face tension and gravity, being capillary and gravity waves in the short and long 
wavelength limits, respectively. As the forcing acceleration is increased, progres- 
sively more complex wave patterns arise, and the inter facial dynamics become 
chaotic [20l [2TJ [19] . Ultimately, large amplitude forcing leads to surface fracture 
and the ejection of droplets from the fluid bath. A recent study on the breaking 
of Faraday waves in a vertically shaken bath has been performed in both the 
capillary and gravity wave limits by Puthenveettil et al. [22. Goodridge et al. 
[23] studied the drop ejection threshold of capillary waves in a glycerine- water 
solution for frequencies up to 100 Hz. 

Faraday [M] reported that such parametric waves can also be emitted by a 
vertical plate plunged into a liquid bath and oscillating horizontally: along both 
sides of the plate, waves aligned perpendicular to the plate oscillate at half the 
forcing frequency. These so-called cross-waves, or edge-induced Faraday waves, 
are typically produced by a wave-maker, and have received considerable atten- 
tion [24[ [25l |26l H?! . Hsieh [29 examined theoretically wave generation in 
a vibrating circular elastic vessel, specifically the axisymmetric capillary waves 
and circumferential crispations that appear in an inviscid fluid subject to radial 
wall displacement. These studies have demonstrated that the excitation of the 
cross waves is due to a parametric resonance. The complexity of this problem 
lies in the nonlinear interactions between the motion of the oscillating rim and 
the resulting wave field. 

Droplets ejected on the liquid surface by breaking Faraday waves may bounce, 
skid and roll before coalescing. A number of recent studies have examined 
droplets bouncing on a vertically vibrating liquid bath below the Faraday thresh- 
old [30^ [31] [32]. The air film between the drop and the liquid surface is squeezed 
and regenerated at each successive bounce, its sustenance precluding coalescence 
and enabling droplet levitation. A similar effect arises on a soap film, a system 
more readily characterized theoretically [35 . The bouncing periodicity depends 
on the size of the drop and the vertical forcing acceleration of the bath [33l [34] . 
Couder et al. [30 have shown that, through the waves emitted at previous 
bounces, some droplets can walk horizontally across the liquid surface. Several 
factors are needed to sustain a so-called "walker" [36]. First, the drop must 
bounce at half the forcing frequency, so that it resonates with the resulting 
Faraday wave field. Second, the bath must be close to the Faraday instability 
threshold so that Faraday waves of large amplitude and spatial extent can be 
excited by the drop impacts. The droplet bounces on the slope of the wave 
emitted at the previous bounce and so receives an impulsive force in a specific 
direction, along which it walks with a constant speed. Such walkers have both 
wave and particle components, and have been shown to exhibit quantum-like 
dynamical behaviour previously thought to be peculiar to the microscopic realm 
[371 EH [Ml [10] • Might such modern physics arise in our ancient bowls? 

The manuscript is divided into two main parts. In §2, we examine the 
acoustics of the Tibetan singing bowls. We begin in §2.1 by reviewing the 
related dynamics and theoretical description of the wine glass [1]. In §2.2, our 
bowls are presented and their deformation spectra analyzed. Then, by adapting 
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the theoretical description of the vibrating wine glass, we infer the Young's 
modulus of the alloy comprising our bowls. In §3, we examine the dynamics of 
flows generated within liquid- filled vibrating bowls. A review of Faraday waves 
and droplet ejection on a vertically shaken bath is presented in §3.1. In §3.2, 
our experimental technique is detailed. In §3.3, we analyze the surface waves 
generated on the liquid bath, and their relation to Faraday waves. In §3.4, we 
examine the limit of large amplitude forcing, in which droplets are ejected by 
breaking Faraday waves. Comparisons are made with experiments performed 
on a vertically shaken liquid bath. Droplet levitation is considered in §3.5, 
where particular attention is given to developing criteria for droplet bouncing 
and exploring the possibility of walking droplets. A summary of our results is 
presented in §4. 

2 Acoustics 
2.1 Background 

Both the wine glass and the Tibetan bowl can be excited by either tapping or 
rubbing its rim. We denote by (n,m) the vibrational mode with n complete 
nodal meridians and m nodal parallels. Tapping excites a number of vibra- 
tional modes [10 , while rubbing excites primarily the (2,0) fundamental mode 
(see Fig. ^p)- Entirely fiexural motion implies radial and tangential displace- 
ments proportional to nsinn^ and cos nO^ respectively, being the azimuthal 
coordinate [41 . For the (2,0) mode the maximum tangential motion is neces- 
sarily half the maximum normal motion. 

A leather mallet can excite bowl vibrations via a stick-slip process, as does 
a finger moving on a wineglass. The moving mallet forces the rim to follow the 
mallet during the stick phase; during the slip phase, the bowl rim relaxes to its 
equilibrium position. This rubbing results in a sound composed of a fundamental 
frequency plus a number of harmonics. While the mallet is in contact with the 
bowl, one of the nodes follows the point of contact [TT], imparting angular 
momentum to the bound liquid. 

To simplify the acoustic analysis, one can approximate the glass or bowl by 
a cylindrical shell with a rigid base and an open top (Fig. [2|l). The system 
can then be described in terms of 7 physical variables, the radius height 
i^o, thickness a. Young's modulus Y and density ps of the cylindrical shell, and 
the frequency / and amplitude A of its oscillating rim. The system can thus 
be described in terms of 4 independent dimensionless groups, which we take 
to be R/H^ A/a, A/R and a Cauchy number Ca = psf'^A'^/Y that indicates 
the relative magnitudes of the inertial and the elastic forces experienced by the 
vibrating rim. 

The sound is emitted by bending waves that deform the rim transversely as 
they propagate. The speed Vb of bending waves on a two-dimensional plate of 
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Figure 2: a) Schematic illustration of a vibrating vessel filled with liquid. Rele- 
vant parameters are the height Hq and the radius R of the vessel, the thickness 
of the rim a, the liquid level the frequency / and the amplitude A of the 
oscillating rim. b) Top view of a vibrating vessel in its fundamental mode (2,0), 
characterized by its 4 nodes and 4 antinodes. 



thickness a is given by (42] 

where Vl is the longitudinal wave speed in the solid (approximately 5200 m/s 
in glass). In order for the bending wave to traverse the perimeter in an integer 
multiple of the period, we require 

1 2TrR 

1 " IT- 

Thus, since V5 ~ vT^? we have / cx a/R^: the frequency increases with rim 
thickness, but decreases with radius. 

A more complete theoretical analysis of the wine glass acoustics [1] can 
be applied to our Tibetan bowls. An ideal cylinder fixed at the bottom is 
considered (Fig. [2^), its wall vibrating with largest amplitude at its free edge or 
rim. The rim's kinetic energy and elastic energy of bending in the mode (2,0) 
are calculated by assuming that the radial position is proportional to cos26>, 
with being the azimuthal coordinate. By considering conservation of total 
energy (kinetic plus elastic bending energy), an expression for the frequency of 
this mode can be deduced: 



When the bowl is partiahy fihed with hquid to a depth H (Fig.[2|L), the frequency 
decreases. French [1] captured this effect by considering the kinetic energy of 
the hquid near the rim, and so deduced the frequency of the fundamental mode: 
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where pi is the liquid density and a ~ 1.25 is a constant indicating the coupling 
efficiency between the rim and fluid displacements. Similarly, frequencies of 
higher modes can be calculated by considering a radial position proportional to 
cos nO and with m nodal parallels [1 : 
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2.2 Tibetan bowls 

Four different antique bowls of different sizes have been studied (Fig. 1). They 
are referred to as Tibet 1, 2, 3 and 4 and their physical characteristics are 
reported in Table 1. 



Bowl name 


Tibet 1 


Tibet 2 


Tibet 3 


Tibet 4 


Symbol 


• 


o 


M 


♦ 


/(2,o) (empty) (Hz) 


236 


187 


347 


428 


Radius R (cm) 


7.5 


8.9 


6 


5.9 


Thickness a (cm) 


0.34 


0.38 


0.31 


0.37 


Mass (kg) 


0.690 


0.814 


0.306 


0.312 


Volume (cm^) 


76 


97 


35 


33 


Density (kg/m^) 


9079 


8366 


8754 


9372 



Table 1: Physical properties of the four Tibetan bowls used in our study. 



When a bowl is struck or rubbed, the sound emitted by the resulting bowl 
vibrations is recorded with a microphone and a fast Fourier transform per- 
formed on the signal. Different peaks clearly appear in the frequency spectrum, 
corresponding to the bowl's different vibrational modes. Fig. [3^ indicates the 
frequency spectrum generated by striking the empty bowl Tibet 4. When the 
bowl is rubbed with a leather wrapped mallet, the lowest mode is excited along 
with its harmonics, an effect known as a mode "lock in" [43 . The frequency 
spectrum of Tibet 4 when rubbed by a leather mallet is presented in Fig. [sJd. 
Due to the bowl asymmetry, two peaks separated by several Hz arise and a 
beating mode is heard. This split is highlighted in a magnification of the first 
peak /(2,o) m the inset of the figure. The deformation shapes are the same with 
both frequencies but there is horizontal angular shift observed to be 7r/4 be- 
tween them for the fundamental modes (2,0) and 7r/2n for other (n, 0) modes. 
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Finally, we note that, owing to the relative squatness of the bowls and the as- 
sociated high energetic penalty of modes with m 7^ 0, only modes (n, 0) were 
excited; henceforth, such modes are denoted simply by n. 
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Figure 3: a) Frequencies excited in the bowl Tibet 4 when struck with a wooden 
mallet. The different peaks correspond to the natural resonant frequencies of the 
bowl, and the associated deformation modes (n, m) are indicated, b) Excited 
frequencies of the bowl Tibet 4 when rubbed with a leather mallet. The first 
peak corresponds to the deformation mode (2,0) and subsequent peaks to the 
harmonics induced by mode lock-in. In the inset, a magnification of the first 
peak provides evidence of its splitting due to the asymmetry of the bowl. 

We can also find a relation between the different mode frequencies. Since the 
speed of the bending wave is proportional to the square root of the frequency, 
V (X v7 since A = v//, we expect A cx Xj^fja- For the mode n, we thus 
have 2i\R = nA^. The frequency of this mode n should thus be proportional to 

/oc^. (6) 

In Fig. [4^, resonant frequencies of the 4 bowls are plotted as a function of their 
corresponding mode n and a power 2 curve is fit onto each curve. In Fig.[4]3, we 
collapse all these curves onto a line by plotting the characteristic speed fR^/a 
as a function of the mode n, thus validating the proposed scaling (|6|. 

It is readily observed that the resonant frequencies decrease when liquid is 
poured into a vessel. In the inset in Fig. [5^, we report the measurements of the 
fundamental frequency of the bowl Tibet 1 as a function of the dimensionless 
liquid height H/Hq. In Fig. [s^i, we report (/o////)^ as a function of (H/Hq)^. 
According to M|, the slope of this curve gives the ratio f In Fig. 5p, we 
present the dependence of (1 + l/n^)(i?^/a)^/^ on (n^ — 1)^ for the deformation 
modes n =2 through 6 of the bowls Tibet 1, 2, 3 and 4. According to ([5|, 
the slope should be equal to jJ-^ — . For each value of the abscissa there are 
4 measurements corresponding to the 4 Tibetan bowls. Data points from the 
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Figure 4: a) Measurements of the resonant frequencies of the different deforma- 
tion modes n for the 4 different Tibetan bowls, b) Characteristic speed fR^ /a 
as a function of deformation mode n for the 4 different bowls. Data are fit by 
a power 2 curve indicating that fnOcn'^a/R'^^ consistent with Eq. [g] 

different Tibetan bowls overlie each other, especially at low n, indicating that 
all bowls have the same value of ratio Y/ pg , and so are likely made of the same 
material. A linear fit gives Y/pg = 8.65 x lO^Pa.m^/kg. 




(H/h/ (n2-1)2 

Figure 5: a) The fundamental frequency of Tibet I's mode (2,0) as a function of 
the water height H/Hq. According to Eq. 0, the slope is equal to f b) The 
dependence of the resonant frequencies / of all the bowls on the deformation 
mode n may be used to infer the ratio Y/pg. According to ([s]), the slope of the 
line shown is 4^f2 from which we infer Y/ps = 8.65 x lO^a.m^/kg. 

We can simply estimate the density for each Tibetan bowl by measuring 
its mass and volume. Masses are measured by a weight scale with an error of 



8 



2 g while bowl volumes are deduced by fluid displacement during submersion. 
The error made on the volume with this method is estimated to be no more 



than 5%. All such density measurements are reported in Table ?? Taking the 
mean value of the densities, equal to 8893 kg/m^, we calculate a Young modulus 
y = 77 ± 6% GPa. This value is in the Young's modulus range of glasses and 
lower than typical brass, copper or bronze alloys for which Y > 100 GPa. 



3 Fluid dynamics 
3.1 Faraday waves 

Consider an inviscid fluid of density p and surface tension a in a horizontal layer 
of uniform depth h in the presence of a gravitational acceleration g. The layer 
is oscillated vertically in a sinusoidal fashion at a forcing frequency /o = a;o/27r, 
amplitude A and acceleration Tg = Ao;^. Above a critical forcing acceleration, 
standing Faraday waves appear on the surface. The associated surface deforma- 
tion a{x^y^t) can be expressed in terms of the container's eigenmodes Sm{x^y) 
as 

a(x, y^t) = ^ am{t)Sm{x, y), (7) 

m 

with am{t) being the oscillating amplitude of the eigenmode m. Benjamin and 
Ursell [16] demonstrate that the coefficients a^(t) satisfy: 

'dm - 2^CO^UJQt)am 0, (8) 

where 



[gkm + ~^mj ^^uh (kmh) (9) 

represents the classic gravity-capillary wave dispersion relation, 

is the dimensionless forcing acceleration, and Bo = is the Bond number. 
km = 27r/Xm is the wave number of the mode m, and Xm is the corresponding 
wavelength. 

In the absence of forcing, F = 0, (|8| describes a simple harmonic oscillation 
with frequency prescribed by ^ corresponding to the free surface vibrations 
of the liquid. When F > 0, ([8| describes a parametric oscillator as the forcing 
term depends on time. The resulting Faraday waves can be either capillary 
or gravity waves according to the magnitude of Bo; specifically, the former and 
latter arise in the respective limits Bo <C 1 and Bo ^ 1. Equation [8] is known as 
the Mathieu equation and cannot be solved analytically since one of the terms 
is time dependent. However, as the forcing is periodic, Floquet theory can be 
applied to show that an inviscid fiuid layer is always unstable to Faraday waves 
with a frequency up that is half the forcing frequency uq = 2ujf [l6|. In the 
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deep water {kh ^ 1), capillary wave {Bo <C 1) limit, we can infer from (|9| a 
Faraday wavelength : 

Xf = {2nY^^{a/pY^Hfo/2)-y^ (11) 

To incorporate the influence of the fluid viscosity, one can add to ([8| a 
phenomenological dissipation term proportional to the velocity [18]: 

dm + 2(3dm +^m(l " 27cosa;o^)am = 0, (12) 

where (3 is the dissipation rate. This dissipation term leads to an acceleration 
threshold for the Faraday instability. Assuming capillary waves in an unbounded 
and infinite depth liquid, the critical acceleration needed to induce parametric 
instability is given by 

Tp^^{p/a)^/^iyoj'f, (13) 

where u is the kinematic viscosity of the fluid |T9j [44] . 

As the forcing amplitude is further increased, the Faraday wave amplitude 
increases progressively until the waves become chaotic. Ultimately, the waves 
break and droplets are ejected from the surface. Since drops will be ejected by 
the breaking of Faraday waves, we expect their diameter to scale as 

d^r^Xr^ {cjlpf"'ff'^ (14) 

in the capillary wave limit. Droplet ejection arises when the destabilizing inertial 
driving force mVg (with m ~ pA|.) exceeds the stabilizing surface tension force 
irXpa. This implies, via ([l4|, a threshold acceleration that scales as: 

r. - lia/py/'f',^' (15) 

The range of validity of these scalings will be investigated in our experimental 
study. 

3.2 Experimental technique 

The experimental setup used for studying the surface waves generated within 
the Tibetan bowls is presented in Fig. |6] A loudspeaker is set in front of the 
bowl, its signal received from a signal function generator then amplified. When 
the applied signal frequency is close to one of the bowl's resonant frequencies, 
it oscillates in the corresponding deformation mode. Conveniently, with this 
method, we can select a single deformation mode. Recall that by striking or 
rubbing the bowls, we excited several modes (Fig.|3|: moreover, rubbing induced 
a rotational motion that followed the mallet. We can now examine the vibration- 
induced flows in a controlled fashion. 

In order to extend the range of natural frequencies, we consider Tibetan 
bowls, wineglasses and soda cans (with their tops removed), whose resonant 
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Loudspeaker 



Figure 6: Experimental setup. The bowl deformation is excited by a loudspeaker 
emitting sound at a frequency corresponding to a particular vibrational mode of 
the bowl. The deformation of the bowl is measured by an accelerometer glued 
to the bowl's outer wall at the height of the liquid surface. 

frequencies span a broad range from 50 to 750 Hz. For each, we can vary 
the resonant frequency by filling it with liquid. The sound emitted by the 
loudspeaker was not powerful enough to induce signifiant oscillations of the 
soda can rim, a problem we eliminated by directly connecting the vibrating 
membrane of the loudspeaker to the rim of the soda can with a rigid rod. 

We measure the acceleration of the rim at an antinode by gluing a lightweight 
accelerometer (PCB-Piezotronics, 352C65) on the bowl's rim at the level of the 
liquid surface. In the following, we characterize the sinusoidal rim oscillation 
by the dimensionless acceleration F defined as the maximal acceleration of the 
rim normalized by the gravitational acceleration F = Acj^/^f. For wineglasses, 
the accelerometer cannot be used since it dramatically alters the resonant fre- 
quency. We thus used a light weight strain gauge, whose effect on the resonance 
frequency is negligible. 

The strain gauge system provides a measurement of the local extension 
length of the rim at an antinode. The length variation of the strain gauge 
is deduced by measuring its electrical resistance with a Wheat stone bridge. To 
deduce the acceleration of the radial rim movement, we deduce a relation be- 
tween the longitudinal extension e and the radial amplitude of the rim A. Then, 
the acceleration F can be readily deduced. The validity of this indirect method 
was tested on Tibet 1. An accelerometer was glued next to a strain gauge at 
an antinode of the bowl Tibet 1. From the strain gauge measurement, we de- 
duced an acceleration that we compared to the direct measurement from the 
accelerometer. These two independent measurements match well for a wide 
range of accelerations. 

French [1 gives a relation for the convex change of length 61 of a curved 
segment of thickness a as a function of the initial mean radius of curvature 
R and the deformed radius of curvature Tc • Si = f ( 7 — where Iq is 
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the initial segment length. Since the deformation is small, we can approximate: 
^ ~ + ^ where x is the radial displacement of the wall pi . The strain gauge 
gives the measurement of the convex longitudinal extension of the rim e = SI/Iq. 
We can thus deduce the maximal radial amplitude of the rim deformation [1 : 

2 



3.3 Surface waves 

Two different liquids were used in our experiment, pure water with density 
= 1000 kg/m^, viscosity v = 1 cSt and surface tension a = 72 mN/m, and 
Dow Corning silicone oil with po = 820 kg/m^, u = 1 cSt and a = 17 A mN/m. 
The fluid depth and resulting natural deformation frequencies of the bowl were 
varied. Specific deformation modes of the bowls were excited acoustically. We 
proceed by reporting the form of the flow induced, specifically, the evolution of 
the free surface with increasing rim forcing. 

In Fig. [7| we present snapshots of the bowl Tibet 1 resonating in its fun- 
damental deformation mode with different T when it is completely filled with 
water. The loudspeaker produces a sinusoidal sound at a frequency /o = 188 Hz 
that corresponds to the mode (2,0) with four associated nodes and antinodes. 
The vibration of the water surface is forced by the horizontal oscillation of the 
rim. When the amplitude of the rim oscillation is small, axisymmetric progres- 
sive capillary waves with frequency commensurate with the excitation frequency 
appear on the liquid surface. Though almost invisible to the naked eye, they can 
be readily detected by appropriate lighting of the liquid surface (Fig. [7^). When 
r is further increased, relatively large amplitude circumferential standing waves 
appear at the water's edge (Fig. [7]3). These standing ripples, aligned perpen- 
dicular to the bowl's edge, are spaced by approximately a Faraday wavelength 



Ai?, as defined in (11). Moreover, their frequency is half that of the axial waves, 
indicating that these waves correspond to classic cross waves or, equivalently, 
edge-induced Faraday waves [M]. More complicated wave modes appear at 
higher excitation amplitude (Fig.]?]:)). At sufficiently high F, the Faraday waves 
break, and water droplets are ejected from the edge of the vessel (Fig. [7]i), 
specifically from the anti-nodes of the oscillating wall. The ejected droplets 
may bounce, slide, and roll on the water surface before eventually coalescing. 

One of our bowls (Tibet 2) resonates readily in both modes (2,0) and (3,0). 
When completely filled with water, the resonant frequencies of its (2,0) and 
(3,0) modes are / = 144 Hz and / = 524 Hz, respectively. In Fig. [sj we observe 
the progression of the surface waves with increasing amplitude for each of these 
modes. Note that for the mode (3,0), since the frequency is higher, the wave- 
lengths are shorter. Moreover, the sound amplitude needed to produce surface 
waves is necessarily higher for mode (3,0) than (2,0). 

The transition from axisymmetric capillary waves to Faraday waves arises at 
a critical acceleration Tf readily measured by the accelerometer. This threshold 
was measured as a function of the forcing frequency, the latter having been 
tuned to excite the fundamental deformation modes of the bowls with different 
liquid levels. We also measured this acceleration threshold for the deformation 
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Figure 7: Evolution of the surface waves in Tibet 1 bowl filled with water and 
excited with a frequency / = 188 Hz corresponding to its fundamental mode 
(2,0). The amplitude of deformation is increasing from left to right: a) F = 1.8, 
b) F = 2.8, c) F = 6.2, d) F = 16.2. 




Figure 8: Evolution of the surface waves in Tibet 2 bowl filled with water, a-d) 
Mode (2,0) is excited with a frequency / = 144 Hz. e-h) Mode (3,0) is excited 
with / = 524 Hz. The amplitude of deformation is increasing from left to right. 
For the deformation mode n = 2, a) F = 1.0, b) F = 1.7, c) F = 5.0, d) 
F = 12.7. For the deformation mode n = 3, e) F = 5.7, f) F = 7.7, g) F = 15.1, 
h) F = 33.8. 



mode (3,0) of the bowl Tibet 2. Higher frequencies were explored with three 
different wineglasses filled to different levels using the strain gauge system. All 
the measurements with silicone oil of viscosity 1 cSt are presented in Fig. 
(lower curve). Consistent with (13), the data suggest a dependence F^ oc j^l 
In Fig. [TO^, we report our measurements of as a function of frequency for 
both silicone oil and distilled water. Each data set is fit by a 5/3 power law. 
Prefactors of 3.5 x 10~^ for water and 1.7 x 10~^ for 1 cSt silicone oil were 
inferred. 
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3.4 Surface fracture 



When the Faraday waves become sufficiently large, they break, leading to droplet 
ejection. A second critical acceleration can thus be measured, T^, above which 
droplets are ejected from the surface. The droplet ejection starts with very 
few droplets ejected, then the ejection rate increases with forcing amplitude. 
Our criterion for onset is that at least two droplets are ejected in a 15 second 
time interval. The dependence of Td on / is presented in Fig. [9] (upper curve) 
for bowls and wine glasses filled with different levels of 1 cSt silicone oil. The 
drop ejection threshold scales as F^^ oc /^/^, which is consistent with the scaling 



presented in (15). 
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Figure 9: Phase diagram indicating the critical acceleration F^? = /S.uj'^ / g above 
which axial capillary waves give way to Faraday waves, and F^^ above which the 
Faraday waves break. The •, o, and ♦ symbols correspond, respectively, to 
the measurements made with the bowls Tibet 1, 2, 3 4 and the ■ to wineglasses 
filled with different levels of 1 cSt silicone oil. The first acceleration threshold 
(lower curve) scales as F^? cx /^/^ and the second (upper curve) as F^^ ex f^^^ . 



In Fig.p^Ojp, our measurements of the dependence of F^^ on / are reported for 
bowls filled with different levels of either 1 cSt silicone oil or water. The two F^^ 



curves collapse when we use the scaling law (15) expected to apply for vertical 
forcing. Specifically, we find 



0.23(<7/p)i/3/'/'. 



(16) 
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The droplet ejection acceleration threshold is thus in accord with measurements 
of Goodridge ^ and Puthenveettil [22 , even though our forcing is horizontal 
rather than vertical. Moreover, our prefactor is consistent with the results of 
both, who reported values between 0.2 and 0.3. 
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Figure 10: Dependence of acceleration thresholds on the frequency / for our 
vessels filled with different levels of pure water (blue symbols) and silicone oil 
of viscosity v = 1 cSt (red symbols). The •, o, Kl and ♦ symbols correspond, 
respectively, to the measurements made with the bowls Tibet 1, 2, 3 4 and the 
■ to wineglasses, a) Each data set of Faraday threshold measurements, F^?, is 



fit by a 5/3 power law consistent with (13). b) The droplet ejection threshold 



F^, is consistent with the scaling (15). 



The diameter of the ejected droplets was measured for several forcing fre- 
quencies. The Tibetan bowls were fully filled with liquid, either 1 cSt silicone 
oil or pure water. Resonant deformation modes were then excited by the loud- 
speaker emitting sinusoidal signals at the appropriate resonant frequency. The 
level of sound was adjusted so that the acceleration of the rim at the antin- 
odes of the bowl were just above the threshold for droplet ejection, F^^. A high 
speed video camera (Phantom) was used to record the ejected droplets, from 
which drop size measurements were taken. We used a liquid-filled glass in or- 
der to extend the frequency range to 720 Hz. Soda cans have very thin walls, 
and very low resonant frequencies. We were thus able to measure droplet sizes 
for the two liquids at a frequency of 98 Hz. Three cumulative distributions of 
ejected droplet sizes are presented on Fig. [TT^ for the bowls Tibet 1, 3 and 4 
resonating in their fundamental deformation modes (2,0). Assuming that these 
distributions are Gaussian, appropriate fits to the cumulative distribution func- 
tions yield the parameters of the Gaussian distribution functions plotted in the 
inset of Fig. pT] 3. The dependencies of the mean droplet size on the forcing fre- 
quency for the two liquids are presented in Fig. [TT]3. Equation (14) adequately 



collapses our data, provided we choose a prefactor of 0.87 . By way of compar- 
ison, Puthenveettil found 0.92 for their experiments with water, and 1.01 with 
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Perfiuoro-compound FC-72 liquid. Donnelly et al. [45] found a prefactor of 0.98 
from his measurements of aerosol water droplets. 




Figure 11: a) Three cumulative distribution functions of ejected droplet sizes 
for the bowls Tibet 1, 3 and 4, fully filled with pure water and resonating in 
their fundamental deformation modes (2,0). The lines are the curves fit by an 
equation of a cumulative distribution function for the Gaussian distribution, b) 
Mean size of the ejected droplets as a function of the frequency of the oscillating 
vessel, either a Tibetan bowl, wineglass or soda can. The two leftmost data 
points are from the measurements made with the soda cans, the two rightmost 
from the wineglass. Error bars are the size of our symbols. The black curve 
indicates a power law dependence with slope -2/3 and prefactor 0.87. The 
dotted line has the same slope, but a prefactor of 1. Inset: The corresponding 
Gaussian functions inferred from the cumulative distribution functions. 



3.5 Bouncing droplets 

With a more viscous fluid (e.g. 10 cSt silicone oil), the waves are less pro- 
nounced, and the fluid is more strongly coupled to the vibrating sidewalls; 
speciflcally, more of the surface oscillates up and down near the wall's anti- 
node. When a droplet of the same liquid is deposited on the surface, it may 
bounce, levitated by the underlying wave fleld. Such sustained levitation was 
not observed in the Tibetan singing bowl with liquid viscosities lower than 10 
cSt. 

In Fig. [l2^, we present a still image of a drop of diameter 0.5 mm bouncing 
on the liquid surface inside the bowl (Tibet 1) resonating at a frequency of 
/o = 188 Hz. The drop has been made by dipping then extracting a syringe 
needle, on the tip of which a capillary bridge forms and breaks, leaving a drop 
that bounces near the oscillating rim. A movie of the bouncing was recorded, 
then vertical slices of each image through the droplet centerline juxtaposed. We 
can thus construct an image illustrating the dynamics of the droplet (Fig, p^). 
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In this case, the drop experience two bounces of shghtly different amphtude 
while the hquid surface (and the rim) oscihates twice. When a smaher droplet 
(of diameter 0.35 mm) is placed on the oscillating surface, the bouncing motion 
can be more complex. In Fig. 12 3, we see that the droplet bounces only once 
during three oscillations of the surface. Fig. 12 i illustrates the corresponding 
trajectory. We note that the bouncing motion of sufficiently small drops can 
become chaotic. 

We sought to sustain walking droplets in our system. Once the horizontal 
forcing amplitude F is just above the Faraday wave threshold F^?, circumferential 
Faraday waves are sustained at the edge of the vessel. Then, when F is increased, 
Faraday waves propagate progressively towards the center of the vessel, their 
amplitude damped by viscosity. Beyond these waves, the liquid surface was 
quiescent unless perturbed by a bouncing droplet, in which case it could sustain 
a field of Faraday waves. This phenomenon was observed in the bowl Tibet 
2 almost fully filled with 10 cSt silicone oil (/o=140 Hz). Just beyond the 
Faraday waves, a bouncing droplet of diameter 500 /im was made such that 
its bouncing frequency corresponded to the Faraday wave frequency, that is, 
half the frequency of the vibrating rim. Such droplets were unable to excite 
sufficiently large Faraday waves to enable them to walk. We note that the usual 
range of walker diameters is between 650 /im and 850 jam [33l |34] : their mass 
is thus 4 times larger than that of our drops. 



4 Conclusion 

We have presented the results of an experimental investigation of the Tibetan 
singing bowl, its acoustics and hydrodynamics. Its acoustical properties are 
similar to those of a wine glass, but its relatively low vibration frequency makes it 
a more efficient generator of edge-induced Faraday waves and droplet generation 
via surface fracture. 

Our observations of the bowl acoustics have been rationalized by adapting 
French's [T| theory of the singing wine glass. This model allowed us to charac- 
terize the bowl acoustics and infer the Young modulus of the alloy constituting 
our antique bowls. The value we found, Y = 77 ± 6% GPa is in the range 
of glass, somewhat smaller than typical brass, copper or bronze alloys. This 
low value of Y and associated low resonant frequency is a critical component 
in the hydrodynamic behavior of singing bowls: bowls with high fundamental 
frequencies are, like the wine glass, relatively inefficient generators of droplets. 

Particular attention has been given to the Faraday waves produced when a 
critical acceleration of the horizontal rim oscillation is exceeded. These have 
been shown to be due to a destabilization of the axial capillary waves similar 
to those observed and studied theoretically |^ . The acceleration thresh- 
old for droplet ejection has also been investigated and rationalized by simple 
scaling arguments. Droplet size was shown to be proportional to the Faraday 
wavelength, and our measurements were consistent with those on a vertically 
shaken liquid surface [22l l23l [45) . 
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Figure 12: a) Still image of a droplet of 0.5 mm diameter bouncing on a surface 
oscillating at 188 Hz, the frequency of the rim. Both droplet and bath are 10 cSt 
Silicon oil. Spatio-temporal diagrams indicate the vertical trajectory of droplets 
of diameter b) 0.5 mm and c) 0.35 mm. Time elapses from left to right and the 
drop's reflection is apparent, d) A single bounce of a droplet of diameter 0.35 
mm is illustrated by an image sequence. Images are spaced by 1 ms. 

We have demonstrated that, following their creation via surface fracture, 
droplets may skip across or roll along the surface of fluid contained within a 
singing bowl. Moreover, careful choice of fluid properties and droplet position 
introduces the possibility of stable bouncing states reminiscent of those on a 
vertically driven free surface [30 . However, stable walking droplets and their 
concomitant quantum behaviour were not observed. Nevertheless, in develop- 
ing hydrodynamic analogues of quantum systems, the edge-forcing examined 
here may be valuable in presenting a lateral gradient in proximity to Faraday 
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threshold. 
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